arXiv:astro-ph/0503216vl 9 Mar 2005 


Mon. Not. R. Astron. Soc. 000. ITI IT!!! ( 2002'! Printed 2 February 2008 


(MN IATeX style file v2.2) 


Effects of External Tidal Field on the Evolution of 
Multi-mass Star Clusters 

K. H. Lee2’2* H. M. Lee^ and H. Sung^ 

^ Astrophysical Research Center for the Structure and Evolution of the Cosmos, Sejong University, Seoul 143-747, Korea 
^ Centre for Astrophysics & Planetary Science, School of Physical Sciences, University of Kent, Canterbury, Kent CT2 7NR, UK 
^Astronomy Program, SEES, Seoul National University, Seoul 151-742, Korea 


2 February 2008 


ABSTRACT 

We present N-body simulations of realistic globular clusters containing initial mass 
function in the galaxy to study effects of tidal field systematically on the proper¬ 
ties of outer parts of globular clusters. Using NBODY6 which takes into account the 
two-body relaxation correctly, we investigate general evolution of globular clusters in 
galactic tidal field. For simplicity, we have employed only spherical components (bulge 
and halo) of the galaxy. Total number of stars in our simulations was about 20,000. 
All simulations were done for several orbital periods in order to understand the devel¬ 
opment of the tidal tails. In our scaled down models, the relaxation time is sufficiently 
short to show the mass segregation effect, but we did not go far enough to see the 
core-collapse, and the fraction of stars lost from the cluster at the end of simulation 
is only about ^ 10%. The radial distribution of extra-tidal stars can be described 
by a power law with a slope around —3. The direction of tidal tails are determined 
by the orbits and locations of the clusters. We find that the length of tidal tails in¬ 
creases towards the apogalacticon and decreases towards the perigalacticon. This is 
an anti-correlation with the strength of the tidal field, caused by the fact that the the 
time-scale for the stars to respond to the potential is similar to the orbital time-scale 
of the cluster. When the length of tidal tails decreases some of the stars in the tidal 
tails are recaptured by the host cluster. From the investigation of velocity anisotropy 
of the model clusters, we find that in the early stages of globular cluster evolution 
the clusters have radial anisotropy in the outermost parts, while clusters are nearly 
isotropic in the cental region. The radial anisotropy decreases with time. 

Key words: stellar dynamics - globular clusters: general - methods: A-body simu¬ 
lations 


1 INTRODUCTION 

Dynamical evolution of globular clusters is induced by a va¬ 
riety of internal and external processes. The main processes 
dominating the evolution of clusters are two-body relaxation 
and tidal shocks. Two-body relaxation makes the velocity 
distribution of stars towards Maxwellian and high-velocity 
stars gain enough energy to escape from the cluster (Spitzer 
& Thunan 1972). The external tidal shocks also accelerate 
the destruction of clusters (see, for example, Gnedin, Lee & 
Ostriker 1999). Clusters experience tidal shocks when they 
pass through the Galactic disk or close to the Galactic bulge. 
The escaped stars may remain in the vicinity of the cluster 
and form tidal tails for a while. Many observational stud¬ 
ies show that the signs of the existence of tidal tails around 


globular clusters (Grillmair et al. 1995; Grillmair et al. 1996; 
Holland, Fahlman & Richer 1997; Lehman & Scholz 1997; 
Leon, Bergond & Vallenari 1999; Leon, Meylan & Combes 
2000; Testa et al. 2000; Odenkirchen et al. 2001, 2003; Sohn 
et al. 2003; Lee et al. 2003, 2004). The features of tidal tails 
are thought to be related to the orbit and location of the 
cluster and the Galactic potential. However the links be¬ 
tween the properties of tidal tails and the cluster dynamics 
have not been fully understood yet. Combes et al. (1999) 
presented extensive N-body simulations of globular clusters, 
and showed the tidal effects quantitatively and geometri¬ 
cally. The main finding of Combes et al. (1999) is that there 
exist two giant tidal tails around the globular cluster along 
cluster’s orbit. They also showed that the density distri¬ 
bution just outside the tidal boundary follows power laws 
to radius, and the extra-tidal stars form clumps which are 
the tracers of the strongest gravitational shocks. However, 
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they used the PM method for the calculation of the grav¬ 
itational potential: this method cannot take into account 
the two-body relaxation effect on the internal evolution of 
the clusters. Yim & Lee (2002) examined the general evolu¬ 
tion of the globular clusters using NBODY6 which properly 
takes into account the effects of two-body relaxation. They 
showed that the density profiles of the clusters appear to 
become somewhat shallower just outside the tidal boundary 
as observed in many clusters and the directions of the tidal 
tails depend on the location in the galaxy as well as the 
cluster orbit. 

Dynamical evolution of many body systems can be stud¬ 
ied by various methods. The most desirable method to study 
dynamical evolution of a globular cluster is a direct N-body 
integration. However, this method becomes almost impos¬ 
sible to use as N becomes large, where N is the number 
of stars, since the computational time increases with . 
The required value of N (of order 10®) is still beyond the 
capability of currently available computers. Fokker-Planck 
method has been commonly used as an alternative to direct 
N-body integration method. The main advantage of Fokker- 
Planck method is that it does not require much computing 
power. However, in this method, the shape of the cluster is 
restricted to sphere (or oblate spheroid for rotating models, 
see Einsel & Spurzem 1999) and the external gravitational 
field can be accounted for only approximately. To investi¬ 
gate tidal tails of globular clusters in the external poten¬ 
tial realistically, we use N-body method in this study, but 
with significantly smaller number of N than realistic globu¬ 
lar clusters. We will focus on the behavior of stars especially 
in the outer parts of the globular clusters. This paper is or¬ 
ganized as follows. In Sect. 2, we briefly describe the method 
and models. And the results of simulations are described in 
Sect. 3. The final section summarizes out major findings. 

2 METHOD AND MODEL 
2.1 NBODY6-I—b 

N-body simulations for this study have been carried out us¬ 
ing NBODY6-I—I- which is a parallelized version of NBODY6 
developed by Aarseth to use parallel super computer 
(Spurzem 1999). We modified the program to include ex¬ 
ternal potential, which is described in §2.2. The code uses 
a direct summation method in computing the gravitational 
potential and adopts Hermite integration scheme which is 
simple to use and highly effective for advancing the single 
particles. The main idea of Hermite integration scheme is to 
employ a fourth-order force polynomial but now the two first 
terms are evaluated by explicit summation over all N parti¬ 
cles, thereby enabling two corrector terms to be formed. The 
code also adopts an Ahmad-Cohen neighbor scheme (Ahmad 
& Cohen 1973) which is based on the idea of separating the 
total force on a particle into two components, 

F — Firr Freg • 

where Firr is the irregular force caused by relatively nearby 
stars and Freg is the regular force from rather distant stars. 
The irregular force varies with time rather rapidly while the 
regular force changes slowly. Since the external force due 
to the host galaxy is expected to vary rather slowly, we re¬ 
garded such force as a regular one. Each part is represented 


by the high-order force polynomial on its own time-step. 
The neighbor force is updated on a time-scale determined 
by the local fluctuations, whereas the smoother regular com¬ 
ponent due to all other members need only be recalculated 
on a longer time-scale, with intermediate values obtained 
by prediction. At each regular time-step a new neighbor list 
is determined using a given neighbor radius for each parti¬ 
cle. It also includes the formation of binaries via three-body 
processes and their interactions with other stars. The main 
force loop in NBODY6 is decomposed in three independent 
DO-loops, advancing all the regularized binaries (1), com¬ 
puting new irregular forces and applying the corrector due 
to them (2), and computing new regular forces, new neigh¬ 
bor lists and applying the regular force corrections (3). The 
entire content of the routines relevant to (2) and (3) have 
been parallelized in NBODY6-I—1-. This is all done fully in 
parallel for the different particles due in this block (Spurzem 
& Baumgardt 2002). Detailed mechanisms of NBODY6 are 
described in Aarseth (1999, 2004) and references therein. 


2.2 The Model Galaxy 

The gravitational potential of the Galaxy is generated by 
stars, gas, and dark matter. There is no single numerical 
method that can describe the evolution of the Galaxy since 
there are many components that react differently. Further, 
we expect that the time-scale for the variation of Galactic 
potential would be much longer than the orbital time-scale of 
the clusters. Therefore the Galactic potential is assumed to 
be independent of time for simplicity. The main components 
of the Galaxy are central bulge, disk, halo, and bar. In the 
present study, the model Galaxy is assumed to be composed 
of two components: halo and bulge, to focus on the effects 
from spherical components only. We have used simple model 
potentials following Lee et al. (1999). The halo component, 
which gives rise to the flat rotation curve at large radii, is 
assumed to have a logarithmic potential, 

^{r)haio = ln(i?c + r^) + const, 

where Rc is the halo core radius and vo is the constant rota¬ 
tion velocity at large r. As for the bulge component, we have 
assumed the Plummer model to allow a steep outward ve¬ 
locity increase from the Galactic centre (Binney & Tremaine 
1987), 

T / \ Gl\^buLge 

Hr)b^lge = - , „ 

Vr2 -I- ri 

where is a parameter that controls the size of bulge. 
The bulges of real galaxies can be better represented by de 
Vaucouleurs’s law, but we have chosen the Plummer 
model for simplicity. Our model have — 230 ~ 240 pc. In 
our numerical calculations all physical quantities are given 
with dimensionless units. The unit of length is chosen to be 
Rsc = 10 kpc. The total mass of the Galaxy within 10 kpc 
is assumed to be Mio = 1.24 x 10^^ Mq (Caldwell & Os- 
triker 1981). The unit of velocity is Vac = (GMio/R)^^^ ~ 
230.7km/s. In these units, the rotation velocity vq at the 
distance of 10 kpc is 0.92 and Mbuige = 0.08, Rc ~ 0.7. 
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Figure 1. Orbits of model clusters. Rp is distance of the peri- 
galacticon, Ra is distance of the apogalacticon, and E is the el- 
lipticity of orbits. 


2.3 Initial Conditions of Globular Clusters and 
Their Orbits 

We have used the King model with Wo = 5 and Wo = 7 
as initial models, where Wo is the scaled central potential, 
which determines the degree of central concentration (King 
1966). The density profiles of all mass components are as¬ 
sumed to be the same at the beginning (i.e. no initial mass 
segregation). We choose the initial mass spectrum to be a 
simple power law 


Table 1. Parameters of the cluster models and orbits. 


Model 

Wo 

Tc 

(pc) 

n 

(pc) 

'^peri 

(kpc) 

E apo 

(kpc) 

E 

Period 

(Myr) 

ml 

5 

2.5 

25 

4 

4 

0 

167 

m2 

5 

2.5 

25 

4 

5.5 

0.16 

242 

m3 

5 

2.5 

25 

4 

8.1 

0.34 

294 

m4 

5 

2.5 

25 

4 

10.1 

0.43 

331 

m5 

5 

1.5 

16 

2.5 

3.8 

0.20 

191 

m6 

5 

1.5 

16 

2.5 

6.2 

0.43 

238 

m7 

5 

1.5 

16 

2.5 

8.1 

0.53 

270 

m8 

7 

0.8 

25 

4 

5.5 

0.16 

242 

m9 

7 

0.8 

25 

4 

8.1 

0.34 

294 

mlO 

7 

0.6 

16 

2.5 

2.5 

0 

105 

mil 

7 

0.6 

16 

2.5 

3.8 

0.20 

191 

ml2 

7 

0.6 

16 

2.5 

8.1 

0.53 

270 


M = 4.2 X 10'^Mq, and the number of stars is Af = 18, 809. 
We have ignored the effects of stellar evolution in order to 
concentrate on the pure dynamical processes and also for 
simplicity. We have carried out a total of 12 simulations with 
different cluster orbits. For models with elongated orbits, we 
defined the ellipticity as, 

^ _ Ra-Rp 
Ra + Rp 

where Ra and Rp are apogalacticon and perigalacticon, re¬ 
spectively. The parameters of models are listed in Table 1, 
and the cluster orbits are shown in Fig. 1. We have removed 
the stars from our simulations when the distance from the 
cluster centre exceeds five times of the tidal radius. The tidal 
radius Vt of the cluster is defined as. 



where Me is the cluster mass, and Ma is the galactic mass 
within radius Ra - Here we have assumed that the tidal radius 
is determined by the tidal force at the apogalacticon. We are 
mostly interested in the development of tidal tails, and thus 
we do not follow the cluster until its destruction. 


2.4 Time-Scale Consideration 

The number of stars of around N ~ 20, 000 used in the 
present study is clearly much smaller than that of realistic 
globular clusters. Different dynamical processes operate on 
different time-scales. The cluster’s orbit depends only on 
the orbital parameters (such as the ellipticity and peri- or 
apo-galacticon distance), and not on N. The orbital time 
of the stars near the tidal boundary is also very close to 
the orbital period of the cluster. However, the dynamical 
evolution such as core-collapse and mass segregation take 
place in two-body relaxation time-scale. Since the relaxation 
time depends on local variables such as stellar density and 
velocity dispersion, a single representative value was defined 
as ‘half-mass relaxation time’ by Spitzer & Hart (1975) as 


dN{m) oc m °‘dm, 

where dN{m) is the number of stars with masses between m 
and m + dm, and we have chosen a = 2.35 which represents 
a Salpeter initial mass function. The range of mass in our 
models is from 1 to 10 M©. The total mass of a cluster is 


^ 1/2 3/2 

trh = 0 - 138 ^ 1 / 2 ^ 1/2 

where rt, is the radius within which the enclosed mass is the 
half of the total mass and m is the mass of the individual 
stars. For clusters with mass function, m is not a well defined 
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Figure 2. Mass functions of model m3 at a few selected epochs. 
Filled circles, MF from the whole stars; filled triangles, inner re¬ 
gion; filled squares, outer region. We can see that mass segregation 
effect increases with time. 

parameter, but we can simply replace it by the mean mass 
although there is no clear reason (see Lee & Goodman 1995 
for more discussion). As we can see from the above equation, 
trh oc . This means that the dynamical evolution of our 
models takes place much faster than real globular clusters 
for a given orbital period. The initial half-mass relaxation 
time of our model is about 3.7 x 10® yr which is much shorter 
than those of real globular clusters. The evaporation time- 
scale is also proportional to the half-mass relaxation time 
(Henon 1961, Lee & Ostriker 1986), and thus we expect 
the development of the tidal tails would be also relatively 
fast. However, since our simulation would last for only up 
to about 10® years, the entire evolution we observe from 
the simulation could be comparable to that of real globular 
clusters. 

The amount of time for stars remaining in the tidal 
tails is expected to be of order of orbital time-scale of the 
cluster (e.g., Lee & Ostriker 1986), which is not related to 
the number of stars. Therefore, the dynamical features in 
the tidal tails of our simulation would not be affected by the 
fact that we used the smaller number of stars than realistic 
clusters. In Table 1, we have listed the radial period which is 
the time from perigalacticon to next perigalacticon of orbits 
in units of years for our models. 

3 RESULTS 

3.1 Mass Segregation 

During the dynamical evolution of the globular cluster, the 
stars in the cluster try to achieve energy equipartition, so 
that all stars have the same kinetic energy. It causes the pro¬ 
cess that higher mass stars slow down and sink towards the 
central region of the cluster, whereas the lower mass stars 


Figure 3. Radial density profile of models m4, m6, m9, and mlO. 
The initial profiles are shown for comparison. Note the clear in¬ 
dication of tidal tails. The average power law gradient is —3.2. 



r/r, 

Figure 4. Fit of power laws to the external profile of model m9. 
Triangle, the profile of massive stars; filled squares, the profile of 
less massive stars; filled circles, the profile of all stars; solid lines 
are the fitting power laws. The profile of less massive stars shows 
a steeper slope. 

have higher velocities than the mean and tend to occupy 
the outer region of the cluster. This process is called mass 
segregation (Spitzer 1987). Many observational studies have 
shown the mass segregation effect in many globular clusters. 
As a result of mass segregation, the envelope of a cluster is 
preferentially populated by low-mass stars, so they escape 
from the cluster more easily than the high-mass ones. Unlike 
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Figure 5. Development of tidal tails of model m5. Each time step is marked as Myr. Two arrows of each panel indicate the directions 
of the galactic centre (long) and its orbit (short). The unit of coordinate is pc. 

our work, Combes et al. (1999) have employed models with tions becomes steeper in the outer part of the cluster as a 

low mass stars being more populated in the outer parts to result of mass segregation, 

mimic the clusters that have reached the mass segregation, 

... j 1 j i • 1 j iu i u j 1 i- Fig. 2 shows the MFs of model m3 at a few selected 

since their models do not include the two-body relaxation ° 

effects. In the present study, we used models which have no ^^e initial models, MFs are the same in all regions, 

initial mass segregation but evolve towards the segregated that the MFs depend 

TTT r on the position in the cluster. We can see that the massive 

state automatically. We can see the process of mass segre- ^ 

gation in the variation of mass functions (MFs) of globular ^^ars are more concentrated to the inner region while the 
clusters with time. The MF can be represented by a power 1®®® ’^''® dominant in the outer region. As a 

law shown in the previous section. The slope of mass func- the MFs of outer region are steeper than those of 

inner region. We can also see that the mass segregation ef¬ 
fect increases with time. Because of preferential migration of 
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high-mass stars towards the centre, the slope of MF from the 
inner region tends to be more flattened with time. However 
the slope of MF from the whole region does not change very 
much until M ~ 0.5Mi, where Mi is initial mass of a clus¬ 
ter (Lee, Fahlman & Richer 1991; Takahashi & Lee 2000). 
Since the main purpose of our simulation is to examine the 
development and evolution of tidal tails, we did not continue 
the simulations towards the significant loss the total mass. 
Thus the MF of the entire cluster does not change signifi¬ 
cantly until the end of our simulations, but the outer parts 
of the clusters become quickly dominated by low mass stars 
since the mass segregation takes place much faster than the 
evaporation time-scale. 


3.2 Radial Density Profile 

Indications of tidal tails around globular clusters have been 
found in several previous studies of globular cluster radial 
density prohles. For a given tidal held, the outer part of a 
globular cluster has more stars than that of a King model. 
As a result, the slope of the radial density prohle shows a 
deviation from the King model around the tidal radius of 
the cluster. We show some of the examples in Fig. 3. The 
initial prohles are shown for comparison. 

The radial distribution of extra-tidal stars can be de¬ 
scribed by a power law E(r) oc M. We derived the average 
value of 7 = —3.2 from the simulations. This value is in well 
accordance with Combes et al. (1999) who deduced a slope 
around 7 = —3 from their N-body simulations. The reason 
for such a steep slope is that the stars in the tidal tails are 
preferentially on radial orbits, and the gravitational poten¬ 
tial is almost Keplerian. It has long been known that the 
density prohle falls like ® in such a case. 

On the other hand, the observed globular clusters tend 
to have shallower slopes. The observed studies of Grillmair 
et al. (1995) and Leon et al. (2000) showed that most of 
the slopes are around —1. Lee et al. (2003) found that the 
slope of Galactic globular cluster M92 is 7 = —1.27 from 
the deep CCD photometry. Combes et al. (1999) explained 
the shallower slope in the observations as the contamination 
of background-foreground populations. Although this must 
be a main source of discrepancy, we could consider other 
reasons. As an example, Lee et al. (2003) showed that the 
slopes of the density profiles for stars in different magnitude 
bins have different values in the outer part. 

Testa et al. (2000) obtained the surface density prohle 
of M92 using plates from the Digitized Second Palomar Sky 
Survey. They htted the extra-tidal profile of the cluster to 
a power law and found 7 = —0.85 ± 0.08. Lee et al. (2003) 
derived 7 = —0.82±0.10 using only the bright stars (18.5 < 
V < 20.5), which correspond with the stars in the study 
by Testa et al. (2000) from surface density profile of M92. 
Whereas, using the faint stars (22 <V < 23.5), they found 
a steeper slope of 7 = —1.31 ± 0.09, which is signihcantly 
different from that of the bright stars. 

We found similar trends in our results of N-body simu¬ 
lations. In Fig. 4, we show an example radial density prohles 
for different masses in model m9. We can see that the prohle 
of massive stars have shallower slope than that of less mas¬ 
sive stars. The prohle of massive stars {2Mq < M < lOM©) 
follows a slope of 7 = —2.4 while that of less massive stars 
{IMq < M < 2Mq) have a much steeper slope of 7 = —3.1. 
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Figure 6. The orbit and location of each step of model m3. Step 
numbers are marked. 


This implies that the difference in slopes of the density pro¬ 
hles in different magnitude bins is a real feature of the clus¬ 
ters. Most of observational studies have used only bright 
stars, because of the sensitivity of the photometry. This 
could be a part of the reason for the shallower density dis¬ 
tribution in observational data. The velocity anisotropy de¬ 
pends on the mass of the stars, and the higher mass stars 
tends to be on more tangential orbits than lower mass stars 
(see, for example, Takahashi & Lee 2000). The stars with less 
radial anisotropy would be populated in shallower prohle. 

Since the tidal tails of globular clusters are preferen¬ 
tially formed by the lowest mass stars (Combes et al. 1999), 
it is difficult to study tidal tails of globular clusters using 
only bright stars. This means that the study of globular 
clusters tidal tail should be performed using stars with the 
lowest possible mass that is detectable. Deep CCD studies 
of other nearby clusters are needed to verify whether the 
mass dependence of the density prohle slope is a common 
phenomenon. 


3.3 Tidal Tails 

The stars unbound from the host cluster do not escape in¬ 
stantaneously, but they slowly drift from the globular clus¬ 
ter and form tidal tails. They can form well-dehned features 
outside the cluster, as observed on the sky (Lee et al. 2003 
and references therein). The development of tidal tails of 
one model (m5) is displayed in Fig. 5. Two arrows indicate 
the direction of galactic centre (long) and cluster orbital di¬ 
rection (short). The maps are projected in the (X,Y) plane 
on which the cluster orbits. We can see the development of 
almost symmetrical tidal tails in two directions. 

The directions of tidal tails are known as a good tracer 
of cluster orbital phase. Combes et al. (1999) presented from 
N-body simulations that there exist two tidal tails around 
the globular cluster along its orbit. In large scale, it seems 
obvious that the tidal tails would be aligned with the clus- 
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Figure 7. Contours of levels of model m3. Each time-step is marked as Myr and step number. The positions of each time step in the 
orbit can be traced by step numbers in Fig. 6. Two arrows of each panel indicate the direction of the galactic centre (long) and its orbit 
(short). Note the change of length of tidal tails. The unit of coordinate is pc. 

ter’s orbit. Recently, Odenkirchen et al. (2001, 2003) de- lar clusters do not cover such a large scale either. In present 

tected clear tidal tails around the Galactic globular cluster study, we focused on a variation of tidal tails of globular 

Palomar 5, extending over an arc of 10° on the sky, cor- clusters along its orbits in a relatively small scale (within 

responding to a projected length of 4 kpc at the distance 5rt). We found that in small scale, the directions of tidal 

of the cluster. They also showed that the tidal streams are tails of globular clusters are not consistent with its orbital 

aligned with the orbit of this cluster. However, not all globu- direction. In Fig. 6 and Fig. 7, we show the positions in the 

lar clusters show such clear tidal tails like Palomar 5. Actu- orbit of one model (m3) and its contours of levels. To see 

ally, Palomar 5 is the first, and only globular cluster whose the direction of tidal tails, we show the 9/9o as a, function of 

tidal tails have been detected at a very significant level of time for model m3 in Fig. 8. Here, 9o is the angle between 

confidence so far. Most of observational studies about globu- the directions of galactic centre and orbital motion, and 9 is 
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Figure 7 - continued 


the angle between the directions of galactic centre and tidal 
tail. We can see that the direction of tidal tails is almost 
median of the directions of galactic centre and cluster orbit 
in all galactic positions. 

From these figures, we also find the relation between 
the length of the tidal tails and the orbital phase of the 
cluster. The cluster which has a circular orbit shows almost 
the same length of tidal tails since the strength of tidal tails 
does not change with time. However in the elliptical orbit, 
the lengths of tidal tails change with time depending on the 
galactocetric distances Rg- The clusters on elongated orbits 
experience variation in tidal field. If the Galactic potential 


behaves like Keplerian, the tidal force is proportional to 
while the logarithmic potential gives rise to the Rq^ depen¬ 
dence. In any case, the amount of the variation in tidal force 
for clusters on elongated orbits would be significant. We may 
simply assume that the tidal tails would be longer when the 
tidal force is stronger. 

However, from Fig. 6 and Fig. 7, we can see that the 
length of the tidal tails increases towards the apocentre 
where the tidal force from the bulge is weakest (near or¬ 
bit number 147 and 173 in Fig. 6) and decreases towards 
the pericentre. Yim & Lee (2002) first noted such behavior, 
and explained that the stars usually do not respond instan- 
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Figure 8. Direction of tidal tails of model m3. 9o is the angle 
between the directions of galactic centre and orbital motion, and 
6 is the angle between the directions of galactic centre and tidal 
tail. The direction of tidal tail is almost median of the direction 
of galactic centre and cluster orbital direction. 



taneously to the external force. Since the stars near the tidal 
boundary have orbital periods similar to the cluster’s orbital 
time around the galaxy, the effect of the strongest tidal force 
appears when the tidal force is weakest. 


3.4 Extra-Tidal Stars 

The stars in the tidal tails can be called as extra-tidal stars. 
We show the fraction of these stars as a function of time for 
all runs in Fig. 9 and Fig. 10. The definition of extra-tidal 
stars is clear only when the globular cluster is embedded 
in a steady potential, which is the case for a circular orbit. 
A common approximation is to assume that a star is un¬ 
bound when it exceeds the conventional tidal radius (Mey- 
lan & Heggie 1997), and N-body models show that this leads 
to consistent results (Giersz & Heggie 1997). However, for 
the elongated orbits, the tidal radius varies along the or¬ 
bit. Keenan (1981) found that the limiting radii of observed 
globular clusters are close to the tidal radius at pericentre. 
In this study, simply the stars which lie outside the tidal 
radius are regarded as extra-tidal stars. For the elongated 
orbits, the tidal radius at pericentre have been used as limit¬ 
ing radius. In calculation of the fraction, the stars eliminated 
from the simulations which went beyond the 5 times of tidal 
radius were also included in extra-tidal stars. 

Generally the fraction increases with time, but the clus¬ 
ters having elongated orbits show clear decrease at some 
points (in models m3, m4, m6, m7, m9 and ml2). In Fig. 10, 
we show the variation of the fraction of extra-tidal stars 
of model m9 which has eccentricity of 0.34. The arrows in 
Fig. 10 denote the pericentre in the orbit of m9. We can see 
that the fraction of extra-tidal stars reaches small peaks at 
the pericentre and then slightly decreases. These points cor¬ 
respond to the time when the length of tidal tails decreases. 
We can see that when the length of tidal tails decreases part 
of the stars included in tidal tails are recaptured by the host 
cluster. 

In addition, from the comparison of m2 and m5 in 
Fig. 9, we can see the effect of strengths of external po¬ 
tential on the rate of evaporation of stars. With the same 
initial cluster model, the one which has closer orbit from the 
galactic centre (m5) shows steeper increase in the fraction 
of extra-tidal stars. Obviously, the clusters in stronger tidal 
fields experience more rapid evaporation of stars. 
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Figure 9. Fraction of extra-tidal stars as a function of time for 
all runs except m9. Note that the models having elliptical orbits 
show clear decrease at some points. 



Time (Myr) 


Figure 10. Fraction of extra-tidal stars as a function of time for 
model m9. The arrows denote the pericentre in the orbit. The 
fraction makes small peaks at the pericentre and then slightly 
decreases. 


3.5 Angular Momentum 

The model clusters used in this study have no angular mo¬ 
mentum at the beginning. However, the probability that the 
globular clusters is formed without initial angular momen¬ 
tum is very small. A seed angular momentum in a proto¬ 
cloud would remain in the globular cluster without dissi¬ 
pation, since the time-scale of cluster formation is much 
shorter compared to globular relaxation times. White & 
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Figure 11. Z component of specific angular momentum of clus- Figure 13. The angular speed of model cluster m3 along the 

ters as a function of time for models ml and m3. The arrows equator at a few selected epochs. The values are normalized by 

denote the pericentre of model m3 having elliptical orbit. the maximum angular speed (orbital angular speed). 
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Figure 12. The angular speed of model clusters. The values 
are normalized by the maximum angular speed (orbital angular 
speed) of each orbit. 3 models are plotted altogether in panel (d). 

Shawl (1987) have measured the projected axial ratio (b/a) 
of 100 globular clusters in the Milky Way, where a and b 
denote semi-major and semi-minor radii, respectively. They 
obtained that mean axial ratio < b/a >= 0.93 ± 0.01. They 
argued that the flattened shape of the clusters can be caused 
by either anisotropy in velocity dispersion or rotation. Since 
the velocity distribution of stars in tidally truncated clusters 
tends to become quickly isotropic (Takahashi & Lee 2000), 
the flattening is likely to be caused by rotation (Combes 


1999). Kinematical data also showed that the flattening 
could be explained by rotation, and that the minor axes are 
nearly coincident with the rotation axes (Meylan & Mayor 
1986). 

Einsel & Spurzem (1999) demonstrated that the influ¬ 
ence of rotation on star cluster evolution is not small. The 
core-collapse time could be accelerated significantly by ini¬ 
tial rotation. Though rotation is a natural initial condition 
from collapse of a star-forming cloud, it is not easy to in¬ 
clude the rotation in most of the existing evolutionary mod¬ 
els of star clusters. Kim et al. (2002) showed that the evolu¬ 
tion of the tidally limited cluster is signihcantly accelerated 
by initial rotation for equal mass models. Recently, from 
the multi-mass Fokker-Planck models, Kim, Lee & Spurzem 
(2004) confirmed that the rotation accelerates core-collapse 
and dissolution of globular clusters even for multi-mass mod¬ 
els. 

The initial angular momentum is expected to disappear 
from the cluster with time by the escape of stars carrying 
some angular momentum. Frenk & Fall (1982) showed that 
the young globular clusters are flatter on average than old 
ones in the LMC. Globular clusters with shorter relaxation 
time are also rounder (Davoust & Prugnel 1990), which sup¬ 
ports the loss of rotation with relaxation. Kim et al. (2002, 
2004) showed that the rotational energy in unit of total en¬ 
ergy decreases monotonically with time in the tidally limited 
rotating clusters. However, the clusters could acquire angu¬ 
lar momentum through tidal interactions with the Galaxy. 
We show the Z component of specific angular momentum 
(angular momentum per unit mass) of clusters as a function 
of time for models ml and m3 in Fig. 11. The arrows in 
Fig. 11 denote the pericentre of model m3 having elliptical 
orbit. The two models have the same initial conditions ex¬ 
cept for their orbits. Fig. 11 shows that the clusters gradually 
acquire angular momentum by tidal torque. The model ml 
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having circular orbit (maintaining closer distance from the 
Galactic centre) has steeper increase in angular momentum. 
In model m3 having elliptical orbit, just after the perigalac- 
tic passages, the cluster angular momentum shows rapid in¬ 
crease. From these results, we can conclude that the Galac¬ 
tic tidal field play an important role in providing angular 
momentum to the clusters. The peaks of angular momen¬ 
tum appear near apocentre where the length of tidal tails is 
maximum. 

The maximum angular speed acquired by tidal torque 
would be the synchronous speed, which means the cluster ro¬ 
tates at the same angular speed as the orbital angular speed. 
We show the angular speed normalized by maximum angular 
speed (orbital angular speed) at the tidal radius of clusters 
as a function of time for three models in Fig. 12. We can see 
that the angular speed increases with time, but our simu¬ 
lations do not last long enough to achieve the synchronous 
rotation. When the results for three difference orbits are 
folded together as shown in Fig. 12(d), the behavior of an¬ 
gular speed evolution appears nearly the same for all. We can 
not find any clear difference between orbits (Fig. 12(d)) pos¬ 
sibly because of short time coverage. In case of the model m3 
(covering the longest time), the angular speed would reach 
the maximum in about 7 times orbital period if it linearly 
increases. 

If a fully relaxed globular cluster has a significant ro¬ 
tation it would be the result of tidal interactions with the 
Galaxy because the initial rotation decrease with time and 
finally disappear (Kim et al. 2002, 2004). Lee et al. (2004) 
suggested that the globular cluster NGG 7492 which has sig¬ 
nificantly flattened shape has been much affected by tidal 
shocks from the relatively flat MF slope and its small mass. 
If the flattened shape of NGC 7492 is caused by its rota¬ 
tion, Galactic tidal field must have given important influ¬ 
ences. When we get the radial velocity distribution, we will 
be able to verify the role of tidal field in rotating of this 
cluster. 

From the radial profiles of the rotational velocity Kim 
et al. (2002, 2004) showed that the rotational velocities of 
initially rotating clusters decreases beyond the half-mass ra¬ 
dius. As an observational result, the rotation of the bright 
Galactic globular cluster w Gentauri is almost solid-body un¬ 
til about 15% of the tidal radius, and then falls off quickly 
(Meylan & Mayor 1986, Merritt, Meylan & Mayor 1997). 
However, we found that the rotating clusters accelerated 
by tidal torque have different shape of radial angular speed 
distribution from initially rotating clusters. In Fig. 13, we 
have shown the distribution of angular speed as a func¬ 
tion of radius for model m3. Although the data are rather 
noisy, the angular speed starts from a differential rotation 
with higher angular speed at larger radii to rigid body ro¬ 
tation. This is due to the fact that the tidal torque is more 
effective at the outer parts, but the inner parts gets the 
angular speed through the relaxation. The rotational en¬ 
ergy for a uniformly rotating stellar system with angular 
speed LO is around MrttC' and maximum value of w ob¬ 
tained by the tidal torque is ~ VrotlRo, where Vrot is the 
velocity of galactic rotation. The ratio of rotational kinetic 
energy to gravitational potential energy can thus become 
TrotlW ~ (■^) I where a is the velocity dispersion 

of the cluster. For a cluster at Rq = 5 kpc, rt = 25 pc, and 
(T = 5 km/s, TrotlW could be about 0.05. Such an amount 
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Figure 14. The radial variation of anisotropy parameter of mod¬ 
els ml and m3, ml has time range of 132 ~ 554 Myr, and m3 
has 133 ~ 1065 Myr. All models show radial anisotropy in the 
outer parts of the clusters and the radial anisotropy decreases 
with time. 


of rotation could produce the dynamical ellipticity greater 
than 0.1 (e.g., Kim et al. 2002). 


3.6 Velocity Anisotropy 

It has long been known that isolated star clusters become 
radially anisotropic after core-collapse due to the scattering 
of stars out of the cluster core (Spitzer 1987). The develop¬ 
ment of velocity anisotropy was investigated in some works 
which used various numerical methods. Most of these cal¬ 
culations confirmed the earlier expectations of generation of 
radial anisotropy in the outer parts of the clusters, although 
the detailed behavior depends on the numerical methods. 
However, actual clusters are embedded in the Galactic tidal 
field which should impose a finite boundary to the cluster. 
The introduction of tidal boundaries gives significant effects 
on the velocity anisotropy because the stars could outflow 
over tidal boundary. Giersz & Heggie (1997) found in N- 
body simulations that star clusters in tidal fields remain 
isotropic during their evolution except for the outer parts 
which become tangentially anisotropic due to the preferen¬ 
tial loss of stars on radial orbits. Takahashi & Lee (2000) 
found in their anisotropic Fokker-Planck simulations that 
the clusters develop stronger radial anisotropy in the outer 
part of the cluster during the early phase and then the devel¬ 
opment of radial anisotropy stops, and finally strong tangen¬ 
tial anisotropy appears as the cluster loses mass. Recently, 
Baumgardt & Makino (2003) showed with N-body simula¬ 
tions that the outer parts of the clusters become rapidly 
tangentially anisotropic. 

Though our simulations do not cover enough time of 
cluster mass loss, it might be interesting to look for signs 
of velocity anisotropies in our runs. Following Binney & 
Tremaine (1986), we define an anisotropy parameter /I as, 
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where Vr and vt are the radial and tangential velocity dis¬ 
persions respectively. They are measured in a non-rotating 
coordinate system in which the cluster centre is always at 
rest at the origin. 

Fig. 14 shows the radial variation of f3 for models m9 
and m3 at various times. All models show radial anisotropy 
in the outer parts of the clusters. The velocity anisotropy is 
expected to be changed with time. We can see that the radial 
anisotropy decreases with time. Extra-tidal stars also show 
radial anisotropy and do not show clear difference from out¬ 
ermost parts of bounded stars. From Fig. 14, we also confirm 
the development of radial anisotropy in the outer parts, and 
that the radial anisotropy decreases with time. Except for 
the outermost parts, clusters are nearly isotropic and there is 
no sign of the radially anisotropic velocity dispersion which 
forms in isolated clusters. Mass segregation could be a rea¬ 
son. Since massive stars are more centrally concentrated, 
they tend to be located at lower energy levels which have 
more isotropic velocity distributions (Baumgardt & Makino 
2003). 


4 SUMMARY AND CONCLUSIONS 

We have carried out a dozen N-body simulations of the glob¬ 
ular clusters within the external tidal field. Our main con¬ 
clusion can be summarized as follows: 

(i) All runs show very quick development of tidal tails 
within a few orbits. 

(ii) Mass segregation effect appears in all runs while the 
initial mass function is still being preserved. Since the rate 
of stellar evaporation depends on the mass of the stars, the 
mass function would change substantially when the amount 
of evaporated stars becomes more than 50% of the initial 
mass. 

(hi) The radial distribution of extra-tidal stars can be 
described by a power law with a slope around —3.2. This 
result is well consistent with previous study (Combes et al. 
1999). In addition, the radial profile of massive stars have 
shallower slope in the outer part of the cluster than that of 
less massive stars as observed in M92 (Lee et al. 2003). This 
could be a part of the reason why the observational data 
have shallower slopes than predicted by simulations. 

(iv) The directions of tidal tails are determined by the 
location of the cluster and its orbits. The length of tidal tails 
increases towards the apocentre and then decreases towards 
pericentre. When the length of tidal tails decreases some 
of the stars in the tails are recaptured by the host cluster. 
In addition, we verified that the external potential plays an 
important role in destruction of a cluster. 

(v) The clusters which have no initial angular momen¬ 
tum gain angular momentum from tidal interactions with 
galaxy. Angular momentum of model clusters gradually in¬ 
creases with time. The model clusters having elongated or¬ 
bits show rapid increases in angular momentum just after 
the perigalactic passages. 

(vi) The rotating clusters accelerated by tidal torque start 
with the rotational angular speed decreasing towards the 
centre but evolve towards the rigid-body rotation. 


(vii) In the early stages of globular cluster evolution, the 
clusters have radial anisotropy in the outermost parts, while 
clusters are nearly isotropic in the central region. And the 
radial anisotropy decreases with time. 

These N-body simulations help to understand the re¬ 
cent observations of extended tidal tails around globular 
clusters. In this study we did not consider tidal shock by 
galactic disk for simplicity. However disk shock also plays 
important role in the evolution of globular clusters (Gnedin 
et al. 1999). Disk component should be contained in sub¬ 
sequent studies to reproduce more realistic Galaxy models. 
Although we used many particles to reproduce real clusters 
as possible, it is still far from the real ones. It is necessary 
to perform N-body simulations with large number of stars 
up to A ~ 10®. 
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